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Abstract. High precision photometry and spectroscopy of low-mass stars reveal a variety of 
properties standard stellar evolution cannot predict. Rotation, an essential ingredient of stellar 
evolution, is a step towards resolving the discrepancy between model predictions and observations. 

The first rotating stellar model, continuously tracing a low-mass star from the pre-main sequence 
onto the horizontal branch, is presented. The predicted luminosity functions of globular clusters and 
surface rotation velocities on the horizontal branch are discussed. 
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INTRODUCTION 



We include rotation into the stellar evolution code FRANEC (IPietrinferni et all 12004). 
The stellar structure equ ations were corrected for the effects of rotation according to 
the prescription given in lKippenhahn & Thomasl ( 1970|) . allowing a ID treatment of the 
stellar structure. 

Additional to chemical element transport via convection and atomic diffusion, ro- 
tational mixing (described only by the shear instability) is introduced. The chemical 
element and angular momentum transport induced by rotational mixing is treated a s 
a diffusive process with a diffusion coefficient given by iDenissenkov & Weissl (|2004|) . 
The change in the total angular momentum is accounted for during th e entire evolutio n 
of the stellar model by disk braking (Pre-Main Sequence, PMS, see IShu et al.lll994l) . 
magnetic-wind braking (Main Sequence, MS, see lChaboyer et al.lll995l) . and mass loss 
(Red-Giant Branch, RGB). 



CALIBRATION 



Realistic initial conditions have to be found and free parameters have to be calibrated 
in order to predict the rotational properties of the stellar model and their effect on the 
evolution. 

The initial angular rotation velocity C0q is chosen to reproduc e the mean rotational pe- 
riod (P\ Myr) of 5 d derived from a sample of young PMS stars (|Rebullll200lLlHerbst et all 
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20021 and lLittlefair et al. 2005|) in the mass range used for calibrating surface-rotation 
rates for the Hyades. To determine the initial rotation rates (Oq for lower metallicities, 
the total angular momentum of stars of the same mass, but different metallicities, is as- 




FIGURE 1. Left: The observ ed (vsint) in the Hyade s from Soderblom et al.1 dl993h as a function of J e ff 
including a ~4.3kms error dSoderblometallll993h as a function of log(r e ff), derived from their (B- 
V) colours with the transformation relationship by lAlonso et al.l dl996l) . The profile of the B e ff functions 
describing the generation efficiency of magnetic fields in the magnetic-wind braking scenario are given 
in the inset, including the point (M W ;B W ) that defines the shape of the Z? e ff function. Isochrones for (vsini) 
are shown using different B eS functions: (0.005;0.500) solid, (0.007;0.500) long dashed, and (0.009;0.508) 
short-dashed line. The isochrone with (0.005;0.500) has the best % 2 minimisation fit to the data and this 
B e g function is chosen for the further work. Right: The rate of angular momentum loss (^) by magnetic- 
wind braking is plotted against the angular surface rotation rate (G>). The rotational properties of the solar 
models with TDB=30Myr, but three different initial angular velocities are shown: (Oq = 4.064 x 10~ 7 s _1 
(Pi M yr= 12 d- light grey line), coq = 9.740 x lO^s" 1 (P 1My r=5 d - grey line), and coo = 4.870 x 10~ 6 s~' 
(PiMyr=l d - black line). Observed 4^ values for the Sun (log = 0) and a theoretically derived upper 
limit of a realistic 2D Weber-Davis solar-wind model dWeber & Davisl 119671) for the MS (dotted) is 
shown. 



sumed to be the same at the start of their life as a PMS star. Disk braking ceases at a time 
^db = 5 Myr after which no significant amount of circum-stellar disks can be observed. 

The free parameters linked to the mixing length and the magnetic-wind braking are 
calibrated comparing models with observed properties of the Sun (e.g. radius, luminos- 
ity, and surface rotation rate) and surface rotation rates observed for the Hyades (see 
left panel of Fig. 1). The influence of changes to the magnetic-wind braking prescrip- 
tion (e.g. the 5 e ff function, see inset on the left panel of Fig. 1) is also demonstrated by 
the solid, long-dashed, and short-dashed lines in the left panel of Fig. 1. The calibra- 
tion of the braking processes is furthermore confirmed by comparing the rate of angular 
momentum loss of our models to the present Sun and an upper limit of a realistic 2D 
Weber-Davis solar-wind model (We ber & Davisl. Il967h (see right panel of Fig. 2). 



RESULTS 



We aim to produce isochrones and Luminosity Functions (LFs) of globular clusters in 
the Galaxy in order to derive the impact of rota tion with our realistic ro tating stellar 
model to improve on prior exploratory work by IVandenberg et ahl (1998). In addition, 




FIGURE 2. Left: The influence of rotation on the CMD of the ZAHB (Z=1(T 3 , Y=0.246) is presented 
for three predicted cases of non-rotating models (solid line), ©o = 0.5 x 10 _6 s _1 (dotted line), and 
O)o= 1.514 x 10 _6 s _1 (dashed line). The ZAHB including the effects of rotation (too = 1.514 x 10~ 6 s~') 
has an increased brightness at (B-V)=0.3 mag (indicated by arrows) of AMy =0.07 mag compared to the 
non-rotating case. Right: The predicted ZAHB surface rotation velocities given as a function of effective 
temperature for a chemical composition on the ZAHB of Z=10~ 3 , Y=0.246. The thick line represents val- 
ues derived with the previously determined Oq. The thin lines mark the upper and lower limit of predicted 
ZAHB rotation v alue s for initially faster or slow er rotat ing stars, respectively. Obs erved H B data in M 15 
dBehr et al.ll2000l and iRecio-Blanco et al.ll2004. M 79 (IRecio-Blanco et all l2004. M 1 3 dPeterson et alJ 
ll995l and lBehr et alj2000h . NGC 2808 dRecio-Blanco et al.ll2004 . and M 80 dRecio-Blanco et all 120041) 
is also included. 



we predict the surface rotation velocity of Zero Age Horizontal Branch (ZAHB) stars. 
The recent increase in observed ZAHB rotatio n rates and our new ro t ating stellar model 
allow improvements to the theoretical work by ISills & Pinsonneaultl (|2000h . All models 
have a chemi cal composition o f Z=10~ 3 and Y=0.246, typical for globular clusters in 
the Galaxy dZinn & WestUl984l) . 

The Colour-Magnitude Diagram (CMD) of ZAHB is determined for varying 0)q 
(non-rotating, 0.5 x 10~ 6 s -1 , and 1.514 x 10~ 6 s -1 ), given in the left panel in Fig. 2. 
Rotation increases the ZAHB brightness, since faster rotating RGB stars develop a more 
massive helium core at the time of the helium flash. 

ZAHB rotation rates are predicted using realistic initial conditions. A relationship on 
the ZAHB between the surface angular velocity, properties of the surface convective 
region, and the amount of angular momentum in the envelope around the helium core 
at the tip of the RGB is determined. It is important to note that the angular momentum 
within only the surface convective region determines the surface rotation rate of the 
ZAHB star. The relationship has been derived with the predicted ZAHB rotation rate of 
a 0.8 M Q model (without mass loss) that has been followed through the helium flash and 
then extended to hotter ZAHB stars. The predicted ZAHB rotation rates for a realistic 
range of C0q determined from the above mentioned sample of young PMS stars are 
compared with observations and given in the right panel of Fig. 2. 

The impact of rotation on LFs was analysed with four stellar models (0.6, 0.7, 0.8, 




FIGURE 3. LFs for initial chemical composition of Z=10~ 3 , Y=0.246 are shown. Left: Two LFs of the 
RGB bump are shown for a 10.5 Gyr old isochrone. The first includes rotation (©o = 1.514 x 10~ 6 s~'; 
solid line) and the second does not (dotted line). Rotation increases the brightness of the RGB bump by 
AMy =0.021 mag. Right: 12 Gyr old LFs are shown. The LF including rotation is plotted as solid lines 
and the non-rotating comparison is plotted as a dotted line. The SGB shows a less steep change in the LF 
when rotation is included, leading to a 1 1 % increase in number of stars at the base of the RGB. 



and 0.85 M Q ) from which two LFs were generated (10.5 Gyr and 12 Gyr). The young 
LF demonstrates that rotation increases the brightness of the RGB bump by 0.021 mag 
(see left panel in Fig. 3). This increase is within the scatter of the brightness difference 
A^'g p between RGB bump and HB for this metallicity dRiello etalll2003h . The old LF 
illustrates how rotation reduces the Sub Giant Branch (SGB) slope and increases the 
number of stars at the base of the RGB by 11 % (right in Fig. 3). The increase in the 
number count at the base of the RGB is still detectable (5 %) when comparing the non- 
rotating 12 Gyr old and rotating 10.5 Gyr old LF, where both cases are in agreement with 
observational constraints. 



SUMMARY 

The rotating models demonstrate the brightening of the ZAHB and RGB bump due to 
rotation. A a>o = 1.514 x 10~ 6 s _1 results in an increased brightness of the ZAHB of 
AMy = 0.07 mag and an increased brightness of the RGB bump of AMy = 0.021 mag 
compared to the non-rotating case. The increase in the helium-core mass and number 
count at the base of the RGB are just two effects of stellar rotation on low-mass stars 
that will be addressed in future work. 

The rotating code and its capability to follow the evolution through the helium flash 
makes it the first stellar evolution code to include rotation during the helium-flash phase 
for low-mass stars. Given realistic initial rotation rates, the code allows us to predict the 
ZAHB velocity distribution. These predicted rotation rates agree better with observa 



tions than rates determined in the most recent theoretical work by ISills & Pinsonneault 



(|2000f) . Our determined ZAHB surface rotation rates explain the apparent "bimodality" 
of the rotation rates at 1 1 000 K (see right panel in Fig. 2) without the need of a bimodal- 
ity in the initial rotation rates and its survival up to the tip of the RGB. 
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